AD-761  149 


INTERACTION  OF  SPHERICAL  SOURCE  FLOW 
AND  AXISYMMETRIC  FREE  JET  WITH  A 
RAREFIED  BACKGROUND 


Prem  K.  Khosla 

Polytechnic  Institute  of  Brooklyn 


Prepared  for: 

Army  Research  Of f i ce - Du r h a m 
May  1972 


DISTRIBUTED  BY: 


National  Technical  Information  Sorvico 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


6HI 9Z0V 


9  9 


AMPA  OHDKR  No.  1442,  Amcndmcnt  a 
nmmm  codc  icao 


na 

* 

ru 

fen 


INTERACTION  OF  SPHERICAL  SOURCE  FLOW 
AND  AXISYMMETRIC  FREE  JET  WITH  A  RAREFIED  BACKGROUND 

by 

PREM  K.  KHOSLA 


>DDC' 

rfppirpnn  arpfr 

JO*  8  I9B  I 

J 


Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U  5  Department  of  Commuct 
Spnnyt.Jd  VA  22131 


POLYTECHNIC  INSTITUTE  OF  BROOKLYN 

DEPARTMENT 

of 

AEROSPACE  ENGINEERING 
and 

APPLIED  MECHANICS 


MAT  1972 


Aamovao  ron  aoaue  rsuaki 
oiamiauTioM  unumitio. 


Vn 

PIBU  BEPBBT  k.  72-25 

P 


Unclassified 

Srnifilv  Classification 


DOCUMENT  CONTROL  DATA  ■  R  &  D 

rS**rufify  c  tms  at  Item!  ion  of  title,  hndy  of  abstract  and  indexing  annotation  must  be  entered  when  the  overattjrcport_  Is^tdassUled) 


,  ORIGINATING  A  C  M  V  I  T  V  ( Corporate  author)  REPORT  SECURITY  CLASSIFICATION 

Polytechnic  Institute  of  Brooklyn  Unclassified 

Dept,  of  Aerospace  Engrg.  &  Applied  Mechanics  |2b.  GROUP 
Route  110,  Farmingdale,  New  York  11735 

3  INTERACTION  OF  SPHERICAL  SOURCE  FLOW  AND  AXISYMMETRIC  FREE  JET 
WITH  A  RAREFIED  BACKGROUND 


4  DESCRIPTIVE  NOTES  (Type  o /  rrporf  end  Inclusive  dates) 

Research  Report _ 

5  au  THOR  IS  I  (First  name,  middle  initial,  last  name) 

Prem  K.  Khosla 


«.  REPORT  OA  TE 

Mav  1972 


fla.  CONTRACT  OR  GRANT  NO. 

DAHC04-69-C-0077 

6.  PROJEC  T  NO. 


7(1.  TOTAL  NO.  OF  PAGES  [7b.  NO,  OF  REFS 


|  9a.  ORIGINATOR'S  REPORT  NUMBER(S) 


PIBAL  Report  No.  72-25 


c.  ARPA  Order  No.  1442,  Amendment  2  [TToTherrEPORt  no(S,  (Any  other  number,  that  may  be  assigned 


d-  Program  Code  9E30 _ _ _ — 

10.  DISTRIBUTION  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


11.  SUPPLEMENTARY  NOTES 


12  SPONSORING  M1LI  T  ARY  ACTIVITY 

U.S.  Army  Research  Of f ice-Durham 
Box  CM,  Duke  Station 
Durham,  North  Carolina  27706 


Ml.  ABSTRACT 


The  interaction  of  a  hypersonic  source  flow  and  an  axisymmetric  free 
jet  with  a  rarefied  background  is  investigated.  Based  on  the  singular- 
nature  of  the  distribution  function,  the  Boltzmann  equation  is  split  into 
two  kinetic  equations.  This  approach  includes  the  basis  for  the  "beam- 
continuum"  method  as. a  special  case.  Application  of  the  simple  moment 
method  leads  to  a  closed  form  solution  for  the  results  arrived  at  by 
Muntz,  Hamel  and  Maguire  for  source  flow.  The  solution  for  the  axisym¬ 
metric  free  jet  also  agrees  with  the  experimental  results  due  to  Brook, 
Hamel  and  Muntz. 


'.,1473 


Unclass i fied 

Srrurm  l  UsMtH  Jti cm 


Unclassi f ied 

Security  Classification 


Free  jet 
Rarefied  gases 
Source  flow 
Hypersonic  flow 


INTERACTION  OF  SPHERICAL  SOURCE  FLOW 


AND  AXISYMMETRIC  FREE  JET  WITH  A  RAREFIED  BACKGROUND 


by 


Prem  K.  Khosla 


This  research  was  supported  by  the  Advanced 
Research  Projects  Agency  of  the  Department 
of  Defense  and  was  monitored  by  the  U.S.  Army 
Research  Of f ice-Durham,  under  Contract 
No.  DAHC04-69-C-0077. 

Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government. 


POLYTECHNIC  INSTITUTE  OF  BROOKLYN 
Department 
of 

Aerospace  Engineering  and  Applied  Mechanics 

May  1972 


PIBAL  Report  No.  72-25 


TABLE  OF  CONTENTS 


Section 


I  Introduction 

II  Basic  Equations 

III  Moment  Method 

IV  Axisymmetric  Jet 
V  Conclusion 

VI  References 


1 

2 

5 

8 

9 

10 


ii 


I.  INTRODUCTION 


The  effect  of  rarefied  backgro-  .id  gas  on  the  hypersonic  source  flow 
has  been  investigated  by  Muntz,  Hamel  and  Maguire^  (MHM)  and  Brook  and 
Kamel2  (BH) .  The  earlier  authors  (MHM)  made  certain  ad  hoc  assumptions 
regarding  the  flow  processes  and  were  able  to  give  a  simple  picture  of 
the  plume  rarefaction.  This  theory  was  modified  by  the  later 
authors  (BH)  ,  who  used  a  kinetic  mode.'  as  represented  by  the  Boltzmann 
equation  and  a  three-fluid  model  comprised  of  jet,  scattered,  and  back¬ 
ground  molecules.  In  their  analysis  the  Boltzmann  equation  was  arbi¬ 
trarily  split  into  three  equations  describing  the  distribution  function 
for  each  type  of  molecule,  and  formal  asymptotic  solutions  were  sought  in 
terms  of  a  rarefaction  parameter. 

In  the  present  report,  we  utilize  the  fact  that  the  distribution 
function  of  the  jet  molecules  is  singular  in  velocity  space.  Subtracting 
this  explicit  singularity  leads  to  the  splitting  of  the  Boltzmann  equation 
into  two  kinetic  equations  whose  asymptotic  solution  for  a  rarefied 
background  is  sought.  A  similar  procedure  by  Grad  has  been  successfully 
applied  to  obtain  the  structure  of  a  strong  shock  wave.  As  one  alterna¬ 
tive  for  carrying  out  a  solution  to  these  equations,  an  approximate 
technique  of  moments,  following  the  well-known  method  of  Mott— Smith,  may 
be  employed.  It  is  shown  that  this  approach  produces  a  conveniently 
simple  closed-form  solution  of  reasonable  accuracy.  Indeed,  it  is  shown 
that  the  previously-cited  au  hoc  theory  of  MHM  is  actually  equivalent  to 
that  due  to  the  moment-met  lod  solution  of  the  present  equations. 

The  theory  is  applied  here  to  two  sample  cases.  that  of  a 
spherical  source  interacting  with  a  rarefied  background  and  that  of  a 
free  axisymmctric  jet  interacting  with  a  rarefied  background. 

It  is  of  further  interest  to  point  out  that  the  present  analysis  pro¬ 
duces  as  a  special  case  the  beam-continuum  equations  derived  by  Rott  and 
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Whittenbury  ,  which  have  been  applied  by  them  to  an  analysis  of  the 
structure  of  a  strong  shock  wave  and  to  the  evaluation  of  sphere  drag  in 

hypersonic  flow.  These  equations  have  recently  been  further  applied  by 

5 

Hot  and  Turcotte  to  the  solution  of  the  flow  over  the  leading-edge  of  a 
sharp  plate.  The  particular  form  of  the  beam-continuum  equations  derivable 
from  the  present  analysis  depends  upon  the  choice  of  the  continuum  part 
of  the  distribution  function.  Some  of  the  alternative  forms  for  the 
distribution  function  are  that  due  to  Mott-Smith  or  Chapman-Enskog .  The 
Mott-Smith  form  produces  a  special  form  of  the  continuum  model  which 
permits  a  closed  form  solution  to  the  problem.  The  Chapman-Enskog  form 
leads  to  the  Navier-Stokes  model  for  the  continuum  portion  and  requires 
a  numerical  solution  of  it.  It  may  be  noted  that  a  systematic  expansion 
based  on  the  Chapman-Enskog  assumption  will  produce  the  present  Mott- 
Smith-based  solution  as  the  leading  approximation. 

Following  Brook  and  Hamel,  we  assume  here  that  the  flow  has  already 
expanded  such  that  the  velocity  has  achieved  its  limiting  value  before 
an  appreciable  effect  of  the  background  is  experienced. 

II.  BASIC  EQUATIONS 

The  Boltzmann  equation  for  hard  sphere  molecular  interactions  can 
be  written  as 

2|.  =  !';r<f'f]'-ff1)gbdbd-:d-  =  J(f,  f) 

=  J(f,  f)  -  «i(r,  ')  f  (1) 

where  f  is  the  dietribution  function,  -  is  the  molecular  velocity  vector, 
g  is  the  relative  speed  of  the  colliding  mole  :ulcs,  v  is  the  collision 
frequency,  b  and  c  arc  two  geometric  variables,  is  the  value  of  f  with 

•  replaced  by  *,  the  velocity  vector  of  the  second  particle.  f'  and  f^ 
are  the  vnlun  of  f  for  v  iocit’cs  of  two  molecules  after  collisions.  We 


write 

f  =  f  j  +  F 


(2) 


where  is  the  distribution  function  for  the  jet  molecules.  Since  the 
jet  has  already  expanued  to  its  limiting  value,  we  may  write  f.  as 

f j  =  (r) 6( — uQ) 5(") M  ')  , 

5  being  the  Dirac  delta  function.  F  corresponds  to  the  continuum  part  of 
the  distribution  and  approaches  the  background  value  at  large  distances 
from  the  source.  Substituting  this  into  Eq.  (1)  leads  to 


vFf.  +  J(f. 


f  . 
D 


) 


Dt  “  ~  VjF  +  2J(f  ,F)  +  J (F ,  F) 

where 

J ( f j  *  F )  =  J f jF^gbdbdcdC 
vF  =  jF^gbdbdedC 

Vj  =  j  f  jgbdbded  r  and  v  =  Vj  +  Vp 


We  nondimensionalize  the  various  quantities  as 


/2RT 


=  Voo  ;  =  R  ;  9/vr^=  ? 


oo 


c,  ~ ) / (2RT,  )  2  =  z,  | 


F  = 


nb 


oo 


( 2RT, 
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i  -  b  -  - 

— 1  =  n  .  ;  —  =  b  ;  —  =  . 

n.*  3  “  ^ 

J 


3 


Thus  Che  Eqs.  (3)  become 


Df  . 
Dt 


v  2  *  * 

V-  n-  r  n  n 


r  41 

=  — °p  !  re  .  —  — •  _  __  1 

(2RV**  ~  +  J<fj'£J>J 


2 £  _  I^oo 

Dt  7T 

>(  2  RT. 


Vroo°2n*r*  r _ nb 


Let 


R  = 
P 


b 

„  2 
nVr  o  n*r* 
op 

( 2RTfa ) ^ 


nwr*  r 

T^r": 


F  +  — r*  +  2j(f 

n 


rF»]  • 


(4) 


;h;re  n; is  the  somc — — -  -  -  —  -iu,  „bTO  and  Tb  are 

the  number  density  and  temperature  of  the  background  gas  and  7  is  a 

representative  relative  collision  velocity,  m  these  equations^  /„* 

and  Rp  are  the  two  parameters,  the  latter  being  proportional  t„  /“se 

nudsen  number.  Since  we  have  assumed  that  the  flow  has  alreadv 

xow  nas  already  expanded 

o  its  final  hypersonic  value,  we  rescale  the  variables  as 


Rn  f •  =  f* 
P  3  3 


R/RP  -  R 


n.R2  =  n 
3  P  j 


Now,  taking  the  limit  n,  /n*-.0  and  r 

boo  P^°°'  we  ^et  to  the  lowest  order, 


Df* 


Dt  2 

R 

P 


DF 


J  =  0 


Dt  -  "  vjP  +  2j(fJ,p) 


(i) 


(ii)  . 


(5) 


Eq-  <5'-li>  Pr°VideS  "  Cl*"  °£  ll-it  to  perfect  vacuum  expan- 

sron.  P.0  can  be  one  of  the  solutions  only  if  the  background  distribution 

function  is  identically  aero.  Even  a  small  value  for  ffc  leads  to  a  finite 
value  of  F. 

in  the  next  section,  a  moment  method  is  presented  to  solve  Eq.  (5). 
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III.  MOfSJT  METHOD 


By  assuming  F  to  be  a  shifted  Maxwellian.  Grad 'was  able  to  predict 
the  structure  of  a  strong  shock  wave.  In  the  present  case,  we  assume 
the  background  distribution  function  F  to  contain  a  variable  numb -  r 
density  and  to  be  of  the  nondimensional  form, 

N. 


F  «= 


jn 


-c 


(6) 


Substituting  this  in  Eq.  tS-ii)  and  multiplying  by  %  and  integrating  over 
the  velocity  space,  we  get 


ijiL(R2V  .  fk  ,  r)d- 


(7) 


R  dR 


The  right-hand  side  of  Eq.  (7)  is  easily  evaluated  by  realizing  the  fact 
that  collisions  must  conserve  momentum  for  the  gas  represented  by  (f^+F), 
So  from  Eq.  (4),  we  get 


2<5J(f*.F)>  -  <?"^Ff^: 


1 7a ) 


fv  f*  term  is  small  as  compared  to  terms  in  Eqs.  (3-i)  or  (4-i)  but  is  of 
1  F  j 

the  same  order  as  other  terms  in  Eqs.  (3-ii)  and  (4-ii) .]  Thus,  Eqs.  (5), 
(7)  and  (7a)  lead  to 


~2  h 

V  N.R  =  const.  =  (j/8  t) 
oo  j 


(8) 


_  _  t  /T  .  T  being  the  stagnation  temperature 

T  stag  b  stag 


and 


,  h  -9  2N> 

~  -mr  (R2N,  ) - r 

FT  dR  D  R 


or 


dN, 


dR 


=  NiVco 


(9) 
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where 


♦ 

T 


U«V 


oo 


-li 


-V* 


CD 


1  +  2V 


CD 


2V 


OO 


erftv^) 


1 


(10) 


Eqs.  (8)  and  (9)  lead  to 

xft.  JL  . 

Nb  dR  R2  F 


the  solution  of  which,  satisfying  the  boundary  conditions  at 
cor.es  out  to  be 

\o  -  exp(-  J/R; 


infinity, 

(11) 


Eq.  (9)  and  the  solution  (li)  are  the  sane  as  obtained  by  HIW.  Use 
of  a  shifted  chapnan-Enskog  distribution  function  for  F  Usds  to  the  un¬ 
known  beam-continuum  equations  of  Rott  and  Khittenbury  which  could  provide 

a  better  approximation  to  the  value  of 

in  order  to  find  the  structure  of  the  distribution  function,  one  has 
to  integrate  Eq.  (5-»i).  In  its  present  form  it  is  very  difficult  to 
solve  Eq.  (5-ii)  without  resort  to  some  sort  of  approximations.  Since 
the  collisional  invariant  moments  of  J(tyF)  have  to  satisfy  relation  (7a) 
we  will  replace  Jltyh  in  (5-ii)  by  ^  <vf  is  further  assumed  to  te 
a  known  function  of  velocity).  Actually  a  better  approximation  could  be 
to  evalua  this  term  by  using  the  solution  of  the  moment  method  given  in 
the  last  section.  Also,  in  order  to  avord  unpleasant  singularities  In 
J(T_,F),  we  will  replace  the  Ansatz  by  the  following: 

n* (r)  -( *-u  )2/2RT.. 

,  .  V  1  ,  e  °  11  M-I'CIT 

(2R-TU)^ 


where  Tn  is  the  constant  value  of  the  temperature  at  which  the  freeze 
tabes  place.  The  splitting  of  m  Boltzmann  equation  still  holds  nrd 
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variation  of  F  i*  still  governed  by  Eq.  (5-ii).  Now,  "j  can  be  written 


as 


_  (2RT  ),j 

‘•j  '  VR>  TT^~ 


oo 


,  -  e  2  1+2  c2 

~  e  +  -575 -  erf  ( C  ) 

,/tt  2C1  1 


NjV  j 


(12) 


where 


2rt, 

C1  =  (2RT“^  (;“Vco 


b  ^ 


c2=  q2+  C,2+  ^ 


Following  Brook  and  Hamel,  Eq.  (5)  can  be  integrated  along  the  charac¬ 
teristic  line  and  we  get 


F  =  f,e 
b 


V.| 

2Rc 


_9 _ 

sin  9 


★  ★ 


N.  e. 

D 


— H - ftan1  ^ - (f-fl)] 

2Rcsin9l  Rsin  9  J 

+  -~r2‘~  "jl  2  J  f~ - r - rT~ - «*p» 

2R  csin  9  1  1+ ( p '/RsinO ) 2 

RcosO  J 


RcosO 

The  last  equation  can  be  put  in  a  simple  form  if  we  write 


P  -  Rsin?tanZ. 


Thus, 


F  = 


v-i  a 

f,  exp  -  — t - : — -] 

b  2Rc  sinS 


*  f* 

3:1  i 

2Rcsin9 


v 
2Rc 


i- 


X_  £ _  „/2 

n  sins  .  /  - v^Z/2Rcsin 9 


j  Nbe 

J'0 


<JZ 


Integration  over  velocity  space  leads  to  a  linear  integral  equation  for 

N^.  An  iterative  solution  is  immediately  suggested  by  the  nature  of  the 

equation.  However,  simple  estimates  can  be  obtained  by  assuming  N  to  be 
•  b 

independent  of  Z.  For  certain  types  of  collisions  this  can  be  a  fairly 
good  approximation  and  has  been  used  by  Kogan^  and  others.  Thus,  we  get 


fb3xP<-  r^r 

2Rc 


sir  9 


-  5  '*[■ 


exp  (- 


2Rc 


-4—t) 
sin  0 


(13) 
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If  the  last  term  in  Eq.  (13)  is  interpreted  as  due  to  scattered  mole¬ 
cules,  we  can  easily  see  that  their  contribution  is  positive  and  of  the 
order  of  (v*/v*).  Apart  from  velocity  dependent  collision  frequency, 
the  present  result  differs  from  that  of  Brook  and  Hamel  in  the  expression 
for  scattered  molecules.  Integration  of  Eq.  (13)  over  velocity  space 
leads  to  the  value  of  N^. 


IV.  AXISYMMETRIC  JET 

Once  again  we  will  assume  that  the  jet  has  already  expanded  to  its 
final  hypersonic  value  before  the  effect  of  background  molecules  becomes 
important.  So  Eq.  (5)  still  applies  to  the  present  problem.  With  a  spher 
ical  coordinate  system  (r,9),  the  solution  of  Eq .  (5-i)  is 

fj  =  n^  (r,  9)  5  ( ,«-uo)  6  (  ri)  6  (  r)  .  <14) 


The  results  of  studies  on  axisymmetric  jets  by  Robertson  and 

Willis6,  and  Grundy^,  indicate  that  the  number  density  is  essentially 

given  by  th^  inviscid  solution  and  the  temperature  freezes  at  different 

values  along  different  streamlines.  For  number  density  n^  one  can  either 

9 

use  the  density  evaluated  from  the  similarity  solution  by  Thornhill  or 
the  co-relation  with  the  characteristics  solution,  due  to  Ashkenas  and 
Sherman10,  can  be  used.  In  the  present  case  the  latter  value  for  is 

used.  Thus, 


n ^ (r ,  9) 


const . 

”1 - 


r 


2 

cos 


A  straightforward  application 
leads  to 


;  m  -  1.365 


■  CD 


of  Mott-Smith  moment 


(15) 

technique  to  Eq.  (5-ii) 
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The  integration  leads  to 


F  2  -I 

N.  *  A  { 'i)cxp(-  —  COS  rr-) 

D  I  R  ‘C 

In  order  to  satisfy  the  boundary  condition  at  R  xd  ,  we  must  have  A.*l. 
Thus, 

★ 

n  2  -  ; 

N.  =  exp  -  —  cos  rr—  (10) 

b 

On  the  jet  centerline  where  0=0,  (10)  reduces  to  Eq.  (11).  So  far  the 
jet  flow  is  hypersonic)  the  '  variation  of  freeze  temperature  docs  not 
influence  Eq.  (16)  to  a  first  order.  Also,  this  simple  expression  for 
variation  of  suggests  that  background  molecules  are  moving  radially 
and  the  penetration  distance  into  the  jet  is  different  along  different 
streamlines.  The  latter  result  is  intuitively  reasonable  as  the  atten¬ 
uation  of  the  background  depends  upon  the  availability  of  the  number  of 
jet  molecules.  As  a  consequence,  the  background  penetrates  farther  along 
streamlines  away  from  the  axis  of  the  jet.  This  solution  agrees  with  the 

■f 

experimental  data  obtained  by  Brook,  Hamel  and  Muntz. 


V.  CONCLUSION 

The  interaction  of  a  hypersonic  source  flow  and  an  axisymmctric  free 
jet  with  a  rarefied  background  has  been  investigated.  Based  on  the  sing¬ 
ular  nature  of  the  distribution  function,  the  Boltzmann  equation  is  split 
into  two  kinetic  equations.  This  approach  includes  the  basis  for  the 
"beam-continuum"  method  as  a  special  case.  Application  of  the  simple- 
moment  mechod  leads  to  a  closed  form  solution  for  the  results  <.riivod 


+An  extended  abstract  of  a  paper  presented  by  Brook,  Hamel  and  Muntz  at 
the  8th  Rarefied  Gas  Dynamics  Symposium,  and  kindly  made  available  by 
Dr.  Brook,  mentiors,  "If  the  background  density  is  assumed  to  decay  expo 
nentially  ,  n^-nv.  exp  (-Rp  ( 0 ) /R)  ,  and  if  we  calculate  R  (■'>),  we  can  see 
-V0Q 1  - x‘*--  and  data  both  agree  with  the  approximation 


that  the  theoretical  results, 
Rn(  h)  =Rn(°)  c°s  1  •  15ri . 
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at  by  MUM  for  source  flow.  The  solution  for  the  nxisymmctric 

free  jet  also  agrees  with  the  experimental  results  duo  to  Brook,  Hamel 

and  Muntz. 
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